3.16.38 \(\int \cos ^3(c+d x) (a+b \sin (c+d x))^2 (A+B \sin (c+d x)) \, dx\) [1538]

3.16.38.1 Optimal result
3.16.38.2 Mathematica [A] (verified)
3.16.38.3 Rubi [A] (verified)
3.16.38.4 Maple [A] (verified)
3.16.38.5 Fricas [A] (verification not implemented)
3.16.38.6 Sympy [A] (verification not implemented)
3.16.38.7 Maxima [A] (verification not implemented)
3.16.38.8 Giac [A] (verification not implemented)
3.16.38.9 Mupad [B] (verification not implemented)

3.16.38.1 Optimal result

Integrand size = 31, antiderivative size = 132 \[ \int \cos ^3(c+d x) (a+b \sin (c+d x))^2 (A+B \sin (c+d x)) \, dx=-\frac {\left (a^2-b^2\right ) (A b-a B) (a+b \sin (c+d x))^3}{3 b^4 d}+\frac {\left (2 a A b-3 a^2 B+b^2 B\right ) (a+b \sin (c+d x))^4}{4 b^4 d}-\frac {(A b-3 a B) (a+b \sin (c+d x))^5}{5 b^4 d}-\frac {B (a+b \sin (c+d x))^6}{6 b^4 d} \]

output
-1/3*(a^2-b^2)*(A*b-B*a)*(a+b*sin(d*x+c))^3/b^4/d+1/4*(2*A*a*b-3*B*a^2+B*b 
^2)*(a+b*sin(d*x+c))^4/b^4/d-1/5*(A*b-3*B*a)*(a+b*sin(d*x+c))^5/b^4/d-1/6* 
B*(a+b*sin(d*x+c))^6/b^4/d
 
3.16.38.2 Mathematica [A] (verified)

Time = 0.17 (sec) , antiderivative size = 111, normalized size of antiderivative = 0.84 \[ \int \cos ^3(c+d x) (a+b \sin (c+d x))^2 (A+B \sin (c+d x)) \, dx=\frac {(a+b \sin (c+d x))^3 \left (-2 a^2 A b+20 A b^3+a^3 B-5 a b^2 B+3 b \left (2 a A b-a^2 B+5 b^2 B\right ) \sin (c+d x)-6 b^2 (2 A b-a B) \sin ^2(c+d x)-10 b^3 B \sin ^3(c+d x)\right )}{60 b^4 d} \]

input
Integrate[Cos[c + d*x]^3*(a + b*Sin[c + d*x])^2*(A + B*Sin[c + d*x]),x]
 
output
((a + b*Sin[c + d*x])^3*(-2*a^2*A*b + 20*A*b^3 + a^3*B - 5*a*b^2*B + 3*b*( 
2*a*A*b - a^2*B + 5*b^2*B)*Sin[c + d*x] - 6*b^2*(2*A*b - a*B)*Sin[c + d*x] 
^2 - 10*b^3*B*Sin[c + d*x]^3))/(60*b^4*d)
 
3.16.38.3 Rubi [A] (verified)

Time = 0.36 (sec) , antiderivative size = 115, normalized size of antiderivative = 0.87, number of steps used = 6, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.161, Rules used = {3042, 3316, 27, 652, 2009}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \cos ^3(c+d x) (a+b \sin (c+d x))^2 (A+B \sin (c+d x)) \, dx\)

\(\Big \downarrow \) 3042

\(\displaystyle \int \cos (c+d x)^3 (a+b \sin (c+d x))^2 (A+B \sin (c+d x))dx\)

\(\Big \downarrow \) 3316

\(\displaystyle \frac {\int \frac {(a+b \sin (c+d x))^2 (A b+B \sin (c+d x) b) \left (b^2-b^2 \sin ^2(c+d x)\right )}{b}d(b \sin (c+d x))}{b^3 d}\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {\int (a+b \sin (c+d x))^2 (A b+B \sin (c+d x) b) \left (b^2-b^2 \sin ^2(c+d x)\right )d(b \sin (c+d x))}{b^4 d}\)

\(\Big \downarrow \) 652

\(\displaystyle \frac {\int \left (-B (a+b \sin (c+d x))^5+(3 a B-A b) (a+b \sin (c+d x))^4+\left (-3 B a^2+2 A b a+b^2 B\right ) (a+b \sin (c+d x))^3+\left (a^2-b^2\right ) (a B-A b) (a+b \sin (c+d x))^2\right )d(b \sin (c+d x))}{b^4 d}\)

\(\Big \downarrow \) 2009

\(\displaystyle \frac {\frac {1}{4} \left (-3 a^2 B+2 a A b+b^2 B\right ) (a+b \sin (c+d x))^4-\frac {1}{3} \left (a^2-b^2\right ) (A b-a B) (a+b \sin (c+d x))^3-\frac {1}{5} (A b-3 a B) (a+b \sin (c+d x))^5-\frac {1}{6} B (a+b \sin (c+d x))^6}{b^4 d}\)

input
Int[Cos[c + d*x]^3*(a + b*Sin[c + d*x])^2*(A + B*Sin[c + d*x]),x]
 
output
(-1/3*((a^2 - b^2)*(A*b - a*B)*(a + b*Sin[c + d*x])^3) + ((2*a*A*b - 3*a^2 
*B + b^2*B)*(a + b*Sin[c + d*x])^4)/4 - ((A*b - 3*a*B)*(a + b*Sin[c + d*x] 
)^5)/5 - (B*(a + b*Sin[c + d*x])^6)/6)/(b^4*d)
 

3.16.38.3.1 Defintions of rubi rules used

rule 27
Int[(a_)*(Fx_), x_Symbol] :> Simp[a   Int[Fx, x], x] /; FreeQ[a, x] &&  !Ma 
tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
 

rule 652
Int[((d_.) + (e_.)*(x_))^(m_.)*((f_.) + (g_.)*(x_))^(n_.)*((a_) + (c_.)*(x_ 
)^2)^(p_.), x_Symbol] :> Int[ExpandIntegrand[(d + e*x)^m*(f + g*x)^n*(a + c 
*x^2)^p, x], x] /; FreeQ[{a, c, d, e, f, g, m, n}, x] && IGtQ[p, 0]
 

rule 2009
Int[u_, x_Symbol] :> Simp[IntSum[u, x], x] /; SumQ[u]
 

rule 3042
Int[u_, x_Symbol] :> Int[DeactivateTrig[u, x], x] /; FunctionOfTrigOfLinear 
Q[u, x]
 

rule 3316
Int[cos[(e_.) + (f_.)*(x_)]^(p_)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_ 
.)*((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)])^(n_.), x_Symbol] :> Simp[1/(b^p* 
f)   Subst[Int[(a + x)^m*(c + (d/b)*x)^n*(b^2 - x^2)^((p - 1)/2), x], x, b* 
Sin[e + f*x]], x] /; FreeQ[{a, b, c, d, e, f, m, n}, x] && IntegerQ[(p - 1) 
/2] && NeQ[a^2 - b^2, 0]
 
3.16.38.4 Maple [A] (verified)

Time = 0.84 (sec) , antiderivative size = 130, normalized size of antiderivative = 0.98

method result size
derivativedivides \(-\frac {\frac {B \,b^{2} \left (\sin ^{6}\left (d x +c \right )\right )}{6}+\frac {\left (A \,b^{2}+2 B a b \right ) \left (\sin ^{5}\left (d x +c \right )\right )}{5}+\frac {\left (\left (a^{2}-b^{2}\right ) B +2 A a b \right ) \left (\sin ^{4}\left (d x +c \right )\right )}{4}+\frac {\left (-2 B a b +\left (a^{2}-b^{2}\right ) A \right ) \left (\sin ^{3}\left (d x +c \right )\right )}{3}+\frac {\left (-2 A a b -B \,a^{2}\right ) \left (\sin ^{2}\left (d x +c \right )\right )}{2}-A \sin \left (d x +c \right ) a^{2}}{d}\) \(130\)
default \(-\frac {\frac {B \,b^{2} \left (\sin ^{6}\left (d x +c \right )\right )}{6}+\frac {\left (A \,b^{2}+2 B a b \right ) \left (\sin ^{5}\left (d x +c \right )\right )}{5}+\frac {\left (\left (a^{2}-b^{2}\right ) B +2 A a b \right ) \left (\sin ^{4}\left (d x +c \right )\right )}{4}+\frac {\left (-2 B a b +\left (a^{2}-b^{2}\right ) A \right ) \left (\sin ^{3}\left (d x +c \right )\right )}{3}+\frac {\left (-2 A a b -B \,a^{2}\right ) \left (\sin ^{2}\left (d x +c \right )\right )}{2}-A \sin \left (d x +c \right ) a^{2}}{d}\) \(130\)
parallelrisch \(\frac {\left (-240 A a b -120 B \,a^{2}-45 B \,b^{2}\right ) \cos \left (2 d x +2 c \right )+\left (80 A \,a^{2}-20 A \,b^{2}-40 B a b \right ) \sin \left (3 d x +3 c \right )+\left (-60 A a b -30 B \,a^{2}\right ) \cos \left (4 d x +4 c \right )+\left (-12 A \,b^{2}-24 B a b \right ) \sin \left (5 d x +5 c \right )+5 B \cos \left (6 d x +6 c \right ) b^{2}+\left (720 A \,a^{2}+120 A \,b^{2}+240 B a b \right ) \sin \left (d x +c \right )+300 A a b +150 B \,a^{2}+40 B \,b^{2}}{960 d}\) \(164\)
risch \(\frac {3 \sin \left (d x +c \right ) A \,a^{2}}{4 d}+\frac {\sin \left (d x +c \right ) A \,b^{2}}{8 d}+\frac {\sin \left (d x +c \right ) B a b}{4 d}+\frac {\cos \left (6 d x +6 c \right ) B \,b^{2}}{192 d}-\frac {\sin \left (5 d x +5 c \right ) A \,b^{2}}{80 d}-\frac {\sin \left (5 d x +5 c \right ) B a b}{40 d}-\frac {\cos \left (4 d x +4 c \right ) A a b}{16 d}-\frac {\cos \left (4 d x +4 c \right ) B \,a^{2}}{32 d}+\frac {A \,a^{2} \sin \left (3 d x +3 c \right )}{12 d}-\frac {\sin \left (3 d x +3 c \right ) A \,b^{2}}{48 d}-\frac {\sin \left (3 d x +3 c \right ) B a b}{24 d}-\frac {\cos \left (2 d x +2 c \right ) A a b}{4 d}-\frac {\cos \left (2 d x +2 c \right ) B \,a^{2}}{8 d}-\frac {3 \cos \left (2 d x +2 c \right ) B \,b^{2}}{64 d}\) \(240\)
norman \(\frac {\frac {\left (8 A a b +4 B \,a^{2}+4 B \,b^{2}\right ) \left (\tan ^{4}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{d}+\frac {\left (8 A a b +4 B \,a^{2}+4 B \,b^{2}\right ) \left (\tan ^{8}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{d}+\frac {2 \left (11 A \,a^{2}+4 A \,b^{2}+8 B a b \right ) \left (\tan ^{3}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{3 d}+\frac {2 \left (11 A \,a^{2}+4 A \,b^{2}+8 B a b \right ) \left (\tan ^{9}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{3 d}+\frac {4 \left (15 A \,a^{2}+2 A \,b^{2}+4 B a b \right ) \left (\tan ^{5}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{5 d}+\frac {4 \left (15 A \,a^{2}+2 A \,b^{2}+4 B a b \right ) \left (\tan ^{7}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{5 d}+\frac {2 \left (2 A a b +B \,a^{2}\right ) \left (\tan ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{d}+\frac {2 \left (2 A a b +B \,a^{2}\right ) \left (\tan ^{10}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{d}+\frac {4 \left (6 A a b +3 B \,a^{2}-2 B \,b^{2}\right ) \left (\tan ^{6}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{3 d}+\frac {2 A \,a^{2} \tan \left (\frac {d x}{2}+\frac {c}{2}\right )}{d}+\frac {2 A \,a^{2} \left (\tan ^{11}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )}{d}}{\left (1+\tan ^{2}\left (\frac {d x}{2}+\frac {c}{2}\right )\right )^{6}}\) \(346\)

input
int(cos(d*x+c)^3*(a+b*sin(d*x+c))^2*(A+B*sin(d*x+c)),x,method=_RETURNVERBO 
SE)
 
output
-1/d*(1/6*B*b^2*sin(d*x+c)^6+1/5*(A*b^2+2*B*a*b)*sin(d*x+c)^5+1/4*((a^2-b^ 
2)*B+2*A*a*b)*sin(d*x+c)^4+1/3*(-2*B*a*b+(a^2-b^2)*A)*sin(d*x+c)^3+1/2*(-2 
*A*a*b-B*a^2)*sin(d*x+c)^2-A*sin(d*x+c)*a^2)
 
3.16.38.5 Fricas [A] (verification not implemented)

Time = 0.29 (sec) , antiderivative size = 120, normalized size of antiderivative = 0.91 \[ \int \cos ^3(c+d x) (a+b \sin (c+d x))^2 (A+B \sin (c+d x)) \, dx=\frac {10 \, B b^{2} \cos \left (d x + c\right )^{6} - 15 \, {\left (B a^{2} + 2 \, A a b + B b^{2}\right )} \cos \left (d x + c\right )^{4} - 4 \, {\left (3 \, {\left (2 \, B a b + A b^{2}\right )} \cos \left (d x + c\right )^{4} - 10 \, A a^{2} - 4 \, B a b - 2 \, A b^{2} - {\left (5 \, A a^{2} + 2 \, B a b + A b^{2}\right )} \cos \left (d x + c\right )^{2}\right )} \sin \left (d x + c\right )}{60 \, d} \]

input
integrate(cos(d*x+c)^3*(a+b*sin(d*x+c))^2*(A+B*sin(d*x+c)),x, algorithm="f 
ricas")
 
output
1/60*(10*B*b^2*cos(d*x + c)^6 - 15*(B*a^2 + 2*A*a*b + B*b^2)*cos(d*x + c)^ 
4 - 4*(3*(2*B*a*b + A*b^2)*cos(d*x + c)^4 - 10*A*a^2 - 4*B*a*b - 2*A*b^2 - 
 (5*A*a^2 + 2*B*a*b + A*b^2)*cos(d*x + c)^2)*sin(d*x + c))/d
 
3.16.38.6 Sympy [A] (verification not implemented)

Time = 0.35 (sec) , antiderivative size = 228, normalized size of antiderivative = 1.73 \[ \int \cos ^3(c+d x) (a+b \sin (c+d x))^2 (A+B \sin (c+d x)) \, dx=\begin {cases} \frac {2 A a^{2} \sin ^{3}{\left (c + d x \right )}}{3 d} + \frac {A a^{2} \sin {\left (c + d x \right )} \cos ^{2}{\left (c + d x \right )}}{d} - \frac {A a b \cos ^{4}{\left (c + d x \right )}}{2 d} + \frac {2 A b^{2} \sin ^{5}{\left (c + d x \right )}}{15 d} + \frac {A b^{2} \sin ^{3}{\left (c + d x \right )} \cos ^{2}{\left (c + d x \right )}}{3 d} - \frac {B a^{2} \cos ^{4}{\left (c + d x \right )}}{4 d} + \frac {4 B a b \sin ^{5}{\left (c + d x \right )}}{15 d} + \frac {2 B a b \sin ^{3}{\left (c + d x \right )} \cos ^{2}{\left (c + d x \right )}}{3 d} - \frac {B b^{2} \sin ^{2}{\left (c + d x \right )} \cos ^{4}{\left (c + d x \right )}}{4 d} - \frac {B b^{2} \cos ^{6}{\left (c + d x \right )}}{12 d} & \text {for}\: d \neq 0 \\x \left (A + B \sin {\left (c \right )}\right ) \left (a + b \sin {\left (c \right )}\right )^{2} \cos ^{3}{\left (c \right )} & \text {otherwise} \end {cases} \]

input
integrate(cos(d*x+c)**3*(a+b*sin(d*x+c))**2*(A+B*sin(d*x+c)),x)
 
output
Piecewise((2*A*a**2*sin(c + d*x)**3/(3*d) + A*a**2*sin(c + d*x)*cos(c + d* 
x)**2/d - A*a*b*cos(c + d*x)**4/(2*d) + 2*A*b**2*sin(c + d*x)**5/(15*d) + 
A*b**2*sin(c + d*x)**3*cos(c + d*x)**2/(3*d) - B*a**2*cos(c + d*x)**4/(4*d 
) + 4*B*a*b*sin(c + d*x)**5/(15*d) + 2*B*a*b*sin(c + d*x)**3*cos(c + d*x)* 
*2/(3*d) - B*b**2*sin(c + d*x)**2*cos(c + d*x)**4/(4*d) - B*b**2*cos(c + d 
*x)**6/(12*d), Ne(d, 0)), (x*(A + B*sin(c))*(a + b*sin(c))**2*cos(c)**3, T 
rue))
 
3.16.38.7 Maxima [A] (verification not implemented)

Time = 0.22 (sec) , antiderivative size = 128, normalized size of antiderivative = 0.97 \[ \int \cos ^3(c+d x) (a+b \sin (c+d x))^2 (A+B \sin (c+d x)) \, dx=-\frac {10 \, B b^{2} \sin \left (d x + c\right )^{6} + 12 \, {\left (2 \, B a b + A b^{2}\right )} \sin \left (d x + c\right )^{5} + 15 \, {\left (B a^{2} + 2 \, A a b - B b^{2}\right )} \sin \left (d x + c\right )^{4} - 60 \, A a^{2} \sin \left (d x + c\right ) + 20 \, {\left (A a^{2} - 2 \, B a b - A b^{2}\right )} \sin \left (d x + c\right )^{3} - 30 \, {\left (B a^{2} + 2 \, A a b\right )} \sin \left (d x + c\right )^{2}}{60 \, d} \]

input
integrate(cos(d*x+c)^3*(a+b*sin(d*x+c))^2*(A+B*sin(d*x+c)),x, algorithm="m 
axima")
 
output
-1/60*(10*B*b^2*sin(d*x + c)^6 + 12*(2*B*a*b + A*b^2)*sin(d*x + c)^5 + 15* 
(B*a^2 + 2*A*a*b - B*b^2)*sin(d*x + c)^4 - 60*A*a^2*sin(d*x + c) + 20*(A*a 
^2 - 2*B*a*b - A*b^2)*sin(d*x + c)^3 - 30*(B*a^2 + 2*A*a*b)*sin(d*x + c)^2 
)/d
 
3.16.38.8 Giac [A] (verification not implemented)

Time = 0.53 (sec) , antiderivative size = 168, normalized size of antiderivative = 1.27 \[ \int \cos ^3(c+d x) (a+b \sin (c+d x))^2 (A+B \sin (c+d x)) \, dx=-\frac {10 \, B b^{2} \sin \left (d x + c\right )^{6} + 24 \, B a b \sin \left (d x + c\right )^{5} + 12 \, A b^{2} \sin \left (d x + c\right )^{5} + 15 \, B a^{2} \sin \left (d x + c\right )^{4} + 30 \, A a b \sin \left (d x + c\right )^{4} - 15 \, B b^{2} \sin \left (d x + c\right )^{4} + 20 \, A a^{2} \sin \left (d x + c\right )^{3} - 40 \, B a b \sin \left (d x + c\right )^{3} - 20 \, A b^{2} \sin \left (d x + c\right )^{3} - 30 \, B a^{2} \sin \left (d x + c\right )^{2} - 60 \, A a b \sin \left (d x + c\right )^{2} - 60 \, A a^{2} \sin \left (d x + c\right )}{60 \, d} \]

input
integrate(cos(d*x+c)^3*(a+b*sin(d*x+c))^2*(A+B*sin(d*x+c)),x, algorithm="g 
iac")
 
output
-1/60*(10*B*b^2*sin(d*x + c)^6 + 24*B*a*b*sin(d*x + c)^5 + 12*A*b^2*sin(d* 
x + c)^5 + 15*B*a^2*sin(d*x + c)^4 + 30*A*a*b*sin(d*x + c)^4 - 15*B*b^2*si 
n(d*x + c)^4 + 20*A*a^2*sin(d*x + c)^3 - 40*B*a*b*sin(d*x + c)^3 - 20*A*b^ 
2*sin(d*x + c)^3 - 30*B*a^2*sin(d*x + c)^2 - 60*A*a*b*sin(d*x + c)^2 - 60* 
A*a^2*sin(d*x + c))/d
 
3.16.38.9 Mupad [B] (verification not implemented)

Time = 12.29 (sec) , antiderivative size = 127, normalized size of antiderivative = 0.96 \[ \int \cos ^3(c+d x) (a+b \sin (c+d x))^2 (A+B \sin (c+d x)) \, dx=\frac {{\sin \left (c+d\,x\right )}^2\,\left (\frac {B\,a^2}{2}+A\,b\,a\right )-{\sin \left (c+d\,x\right )}^5\,\left (\frac {A\,b^2}{5}+\frac {2\,B\,a\,b}{5}\right )+{\sin \left (c+d\,x\right )}^3\,\left (-\frac {A\,a^2}{3}+\frac {2\,B\,a\,b}{3}+\frac {A\,b^2}{3}\right )-{\sin \left (c+d\,x\right )}^4\,\left (\frac {B\,a^2}{4}+\frac {A\,a\,b}{2}-\frac {B\,b^2}{4}\right )-\frac {B\,b^2\,{\sin \left (c+d\,x\right )}^6}{6}+A\,a^2\,\sin \left (c+d\,x\right )}{d} \]

input
int(cos(c + d*x)^3*(A + B*sin(c + d*x))*(a + b*sin(c + d*x))^2,x)
 
output
(sin(c + d*x)^2*((B*a^2)/2 + A*a*b) - sin(c + d*x)^5*((A*b^2)/5 + (2*B*a*b 
)/5) + sin(c + d*x)^3*((A*b^2)/3 - (A*a^2)/3 + (2*B*a*b)/3) - sin(c + d*x) 
^4*((B*a^2)/4 - (B*b^2)/4 + (A*a*b)/2) - (B*b^2*sin(c + d*x)^6)/6 + A*a^2* 
sin(c + d*x))/d